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PREFACE 


The  study  reported  here  was  authorized  by  Headquarters,  U.S.  Army  Corps  of 
Engineers  (HQUSACE),  as  part  of  the  Concrete  and  Steel  Structures  Problem  Area 
of  the  Repair,  Evaluation,  Maintenace,  and  Rehabilitation  (REMR)  Research 
Program.  The  work  was  performed  under  Work  Unit  3301,  "Structures  Damage 
Index  to  Determine  Remaining  Life  and  Reliability  of  Metal  Structures."  Mr. 
Frank  Kearney  was  Principal  Investigator  for  this  Work  Unit.  Dr.  Tony  C.  Liu 
(CECW-ED)  was  the  REMR  Technical  Monitor  for  this  work. 

Mr.  Jesse  A.  Pfeiffer,  Jr.  (CERD-C)  is  the  REMR  Coordinator  at  the 
Directorate  of  Research  and  Development,  HQUSACE;  Mr.  John  R.  Mikel  and  Mr. 
James  Crews  (CECW-OM)  and  Dr.  Tony  C.  Liu  (CECW-ED)  serve  as  the  REMR  Overview 
Committee;  Mr.  William  F.  McCleese  (CEWES-SC-A) ,  U.S.  Army  Engineer  Waterways 
Experiment  Station  (WES),  is  the  REMR  Program  Manager.  Mr.  J.  McDonald 
(CEWES-SC-R)  is  the  Problem  Area  Leaner. 

This  study  was  conducted  by  the  Structures  Division  of  JAYCOR  for  the 
Engineering  and  Materials  Division  (EM)  of  the  U.S.  Army  Construction  Engin¬ 
eering  Research  Laboratory  (USACERL),  during  the  period  April  1985  to 
September  1985.  Appreciation  is  expressed  for  the  assistance  of  Dr.  Dawn 
White,  Mr.  Barry  Hare,  and  Dr.  Anthony  Kao  (all  of  USACERL)  and  Ms.  Caroline 
Cummins,  Ms.  Susan  Lanier,  and  Mr.  William  J.  Flathau  (all  of  JAYCOR).  COL 
Carl  0.  Magnell  is  Commander  and  Director  of  USACERL,  and  Dr.  L.  R.  Shaffer  is 
Technical  Director.  Dr.  Robert  Quattrone  is  chief  of  USACERL-EM. 
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ASSESSING  THE  RELIABILITY  OF  STEEL  CIVIL  WORKS  STRUCTURES 


PART  I:  INTRODUCTION 


Background 


1.  A  major  mission  of  Che  U.S.  Army  Corps  of  Engineers  (USACE)  is 
maintaining  the  operational  efficiency  of  various  navigation,  fLood  control, 
and  hydroelectric  power  projects.  The  importance  of  this  mission  increases 
with  time  as  older  projects  deteriorate  and  few  new  projects  are  authorized. 

Of  particular  concern  are  steel  structures.  A  recent  study  (JAYCOR  1985) 
documented  that  USACE  Civil  Works  projects  include  many  structurally  sig¬ 
nificant  steel  features  such  as  bridges,  outlet  gates,  spillway  gates,  lock 
gates  and  valves,  sheet  piles,  hydropower  gates,  and  penstock  liners  (Figure 
1).  The  study  further  estimated  that  each  year  these  steel  structures  experi¬ 
ence  about  44  different  problems  to  varying  degrees. 

2.  Research  has  shown  that  a  relatively  significant  fraction  of  the 
problems  occur  with  steel  sheet  piling.  Deterministic  analyses  of  this 
structural  feature  have  been  conducted  by  the  U.S.  Army  Construction  Engineer¬ 
ing  Research  Laboratory  (USA-CERL)  (Kearney  1986).  That  study  of  sheet  pile 
structures  identified  the  elements  that  could  fail  (Figure  2). 

3.  A  structural  condition  index  would  make  it  easier  to  allocate  public 
funds  where  they  are  most  needed  to  solve  these  problems  and  avert  failures. 
Such  an  index  should  indicate  the  composite  structural  integrity  and  the 
overall  operational  condition  of  a  particular  feature,  and  reflect  the  extent, 
severity,  and  type  of  deterioration.  It  should  also  take  into  account  the 
effects  of  loads  imposed  on  the  structure.  The  reliability  function  can 
provide  such  a  condition  index.  By  quantifying  changes  in  resistance  and 
loading,  this  function  indicates  the  probability  that  a  structure  will  perform 
adequately  under  certain  conditions.  Also,  as  used  in  this  report,  it  pro¬ 
vides  more  information  than  a  deterministic  analysis  because  it  shows  the 
gradual  deter iorat ion  of  condition  over  time. 
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Figure  2.  Structurally  critical  elements  of  steel  sheet  pile  structures 
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PART  1 1 :  METHODOLOGY 


Reliability  as  an  Index  of  Condition 

7.  This  chapter  describes  a  general  procedure  for  estimating  the 
reliability  of  existing  CE  structures.  The  reliability  of  a  structure  is  the 
probability  that,  when  operating  under  stated  environmental  conditions,  it 
will  perform  its  intended  function  adequately  for  a  specified  interval  of  time 
(Kapur  and  Lamberson  1977).  The  goal  of  reliability  analysis  is  to  find  the 
likelihood  that  the  load  exceeds  the  resistance.  Thus,  each  calculation  has 
two  distinct  parts  for  each  mode  of  failure:  stress  (load)  analysis  and 
strength  (resistance)  analysis.  This  reliability  function  provides  all  the 
characteristics  desirable  for  an  index  of  a  structure's  present  condition.  As 
a  mathematical  "probability,"  it  is  a  numerically  ordered  measure,  with  unity 
representing  the  best  condition  (best  performance)  and  zero  representing  the 
worst  condition  (no  performance).  "Stated  environmental  conditions"  expli¬ 
citly  incorporate  the  effects  of  imposed  service  'oadings  that  were  not 
contemplated  during  design.  The  "specified  interval  of  time"  is  necessary  for 
quantifying  any  deterioration  in  strength  capacity  that  a  structure  experi¬ 
ences.  Finally,  the  "intended  function"  makes  it  possible  to  address  manage¬ 
able  operation  and  maintenance  problems. 

8.  Probabilistic  methods  will  be  most  useful  in  evaluating  existing 
structures,  rather  than  in  designing  new  ones,  because  they  are  better  for 
handling  uncertainties.  For  design,  the  Corps  of  Engineers  has  traditionally 
used  deterministic  rather  than  probabilistic  methods.  For  example,  in 
designing  a  lock  chamber  with  sheet  pile  cells  as  walls,  uncertainties  will 
arise  about  operational  loadings  and  material  strengths.  A  prudent  designer 
uses  deterministic,  conservative  values,  possibly  overdesigning  the  feature. 
The  additional  cost  incurred  by  overdesigning  is  usually  only  a  small  fraction 
of  the  total  construction  cost.  On  the  other  hand,  when  evaluating  a  struc¬ 
ture's  need  for  rehabilitation,  the  same  uncertainties  will  arise  as  in 
design,  along  with  others  arising  from  the  evaluator's  inability  to  make 
certain  measurements  and  inspect  certain  parts.  However,  in  rehabilitation, 
the  costs  of  uncertainty  are  much  greater.  Not  only  is  the  actual  work 
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CURRENT 


Figure  3.  Miter  gates:  horizontally  framed  typical  girder  data 


considered  random  variables.  For  them,  a  mean  and  a  standard  deviation  are 
assigned.  Those  two  parameters  are  precisely  defined  in  elementary  probabil¬ 
ity  texts  (Benjamin  and  Cornell  1970).  The  mean  measures  the  central  tendency 
of  the  variable  or  the  value  about  which  scatter  can  be  expected  if  repeated 
observations  of  the  phenomenon  are  made.  The  standard  deviation,  which  has 
the  same  physical  units  as  the  mean,  measures  the  dispersion  of  repeated 
observations  abput  the  mean  value. 

14.  The  engineer  can  draw  upon  extensive  sources  of  information  to  assign 
values  for  these  parameters.  In  some  cases,  a  sample  of  repeated  observations 
having  useful  descriptive  statistics  may  have  been  collected.  Often,  data 
collected  for  similar  problems  may  apply.  An  experienced  engineer  can 
quantify  his  or  her  expert  judgment  by  selecting  the  mean  and  the  standard 
deviation  for  random  variables.  This  step  is  different  from  the  design 
situation  in  which  a  single  conservative  value  is  assigned  to  uncertain 
factors.  However,  in  evaluating  existing  structures,  it  is  important  to 
quantify  what  is  known  about  variables  by  using  means,  and  to  quantify  what  is 
not  known  about  them  by  using  standard  deviations. 

15.  In  the  example  shown  in  Figure  3,  the  geometry  of  the  girder  is  pre¬ 
sumed  to  be  known  with  confidence.  Accordingly,  the  values  for  these  para¬ 
meters  will  be  deterministically  assigned,  as  given  in  Table  1.  These  numbers 
correspond  to  an  example  design  (Granade  1980),  except  for  the  moment  of 
inertia,  I,  which  is  reduced  by  30  percent  to  reflect  hypothetical  corrosion 

of  the  upstream  flange.  The  yield  stress  of  the  structural  steel,  fy,  is 
assumed  to  be  random.  It  is  reasonable  to  expect  that  the  findings  of  a 
general  study  (American  Iron  and  Steel  Institute  1978)  of  this  factor's  varia¬ 


bility  apply  to  the  example  outlined  here.  Accordingly,  the  mean  and  standard 
deviation  of  f  are  assigned  as  indicated  in  Table  1.  Note  that  the  mean 
value  is  greater  than  the  36,000  psi  normally  used  in  design.  The  girder 
loading  is  also  taken  to  be  random.  Judgment  is  used  to  assign  its  mean  value 
to  be  150  percent  of  the  hydrostatic  value,  reflecting  the  effects  of  ice, 
mud,  and  live  loads.  The  uncertu'nty  about  this  factor  is  quantified  in  a 
standard  deviation  of  20  percent  of  the  mean. 


of  these  results  using  Equations  4  and  5.  More  generally,  if  MS  (X,,  X2 ,  ... 
Xn)  is  a  function  of  n  independent  input  variables: 


where  ms  is  the  margin  of  safety  calculated  using  the  mean  value  of  each 
variable: 

ms  =  MS (X  ,  X,  ....  X  )  (10) 

12  n 

anti  MS^  and  0  ate  computed  from  Equations  4  and  5  as  if  X-  were  the 

only  random  variable  and  the  others  were  equal  to  their  mean  values.  For  this 
general  case,  2n  +  1  design-likc  evaluations  of  the  margin  of  safety  are 
required . 

18.  With  these  estimates  of  MS  and  c  ,  a  convenient  assumption  can  be 
made  about  the  form  of  the  MS  distribution  in  order  to  complete  the  estimation 
of  reliability.  Since  the  margin  of  safety  is  the  difference  between  resis¬ 
tance  and  load  which,  in  turn,  often  involves  sums  of  random  variables,  there 
is  some  justification  to  approximate  MS  with  a  normal  distribution  (Benjamin 
and  Cornell  1970).  Since  the  structural  evaluations  of  operational  interest 
do  not  involve  extremely  rare  events,  in  most  cases  the  reliabilities  calcu¬ 
lated  will  not  be  extremely  sensitive  to  the  form  of  this  distribution.  The 
assumption  of  a  normal  distribution  has  tuo  other  advantages:  it  is  familiar 
even  to  those  new  to  probability,  and  standard  tables  and  approximating 
functions  are  available.  Finally,  it  follows  that  the  reliability  is  given 
by: 

R  =  Probability  (MS  >  0)  =  1  -  F  (-  —  |  (11) 

\  °MS  / 

in  which  F()  is  the  standard  normal  cumulative  distribution  function  (Benjamin 
and  Cornell  1970). 
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Reliability,  wh:  indexes  structmal  condition  on  a  zero  to  one  scale, 

follows  from  Equation  11.  Note  that  this  sequence  is  identical  irrespective 
of  the  deterministic  model  selected  in  Step  1.  It  applies  for  a  comprehen¬ 
sive,  finite  element  stress  analysis  as  well  as  for  a  statically  determinate 
equilibrium  equation.  Appendix  A  provides  a  microcomputer  program  to  imple¬ 
ment  this  procedure. 


Exam pie  Calculation 

21.  For  the  example  given  in  Figure  3,  n  =  2  variables  have  beer:  ::n- 
sidered  random,  f  and  w.  Accordingly,  five  deterministic  evaluations  or 
Equation  3  are  required,  as  summarized  in  Table  2.  The  estimates  of  MS  and 
gMS’  cons r i n8  only  f  to  be  random,  are  from  Equations  4  and  5  and  the 
second  and  third  rows  of  Table  2: 


MS 


MS 


MS ,  = 

i 


4271  +  13071 

- - -  -  86  71  ps  i 


MS  -  MS 
+ 


13071  -  4271 


MS. 


=  4400  psi 


Considering  only  w  to  be  random,  the  last  two  rows  of  Table  2  gave: 


MS 


2 


a 


MS 


2 


14856  +  2485 
2 

14856  -  2485 
2 


8671  psi 
6186  psi 


Substitution  of  these  results  and  ms  from  the  first  row  of  Table  2,  into  Eq  8 
leads  to: 


MS 

867  1 


=  (1.000)  (1.000) 


t  rom  which  MS  =  867  1  psi.  In  turn.  Equation  9  becomes: 


/  ^uc  \  y  9 

1  +  18671/  =  (1  *  °-50/')  (1  +  0.713  ) 


°MS  =  8208  psi 
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PART  III:  O'BRIEN  LOCK  AND  DAM 


24.  In  1979,  an  inspection  of  the  O'Brien  Lock,  and  Dam,  located  on  the 
Illinois  waterway  in  Chicago,  revealed  excessive  corrosion  of  sheet  piling.  A 
deterministic  analysis  of  the  problem  (Kearney  1986)  concluded  that  the  struc¬ 
ture  should  have  a  service  Life  of  83  years.  For  this  study,  the  O'Brien  Lock 
and  Dam  has  been  used  as  an  example  of  how  probabilistic  analysis  can  be 
employed  to  calculate  service  life. 

General  Information 


25.  The  O'Brien  Lock  and  Dam,  a  cellular  cofferdam,  presents  a  complex 
problem  because  the  factors  affecting  its  service  life  are  many  and  subtle. 
Cellular  cofferdams  are  double-wall  retaining  structures  constructed  of 
interlocking  steel  sheet  piles.  The  piles  form  adjacent  cells  that  are  filled 
with  soil  or  rock  fragments  (Lacroix,  Esrig,  and  Luscher  1970).  In  analyzing 
cofferdams,  one  needs  to  consider  not  only  the  structure's  individual  charac¬ 
teristics,  but  also  the  interactions  among  soil,  rock,  water,  weather,  steel, 
and  concr°te.  A  dam's  shape  and  location  also  make  it  an  individual  system  of 
unique  characteristics:  cofferdams  may  be  shaped  in  a  circle,  diaphragm,  or 
cloverleaf;  sheet  piling  walls  can  be  constructed  on  dry  land,  in  fresh  water, 
or  in  seawater.  However,  all  cofferdams  have  one  factor  in  common:  any 
interlock  in  a  single  sheet  pile  cell  is  a  possible  failure  location,  and 
unlike  most  conventional  structures,  there  is  little  redundancy  in  the  system. 

26.  The  following  types  of  cofferdam  failures  have  been  observed  and 
reported : 

a.  Sliding  on  the  base. 

b.  Bank  sliding. 

c.  Flooding  of  a  cell. 

d.  Erosion  in  streams — scour  and  loss  of  cell  fill. 

e.  Boiling  of  land  cofferdams. 

f.  Shearing  of  silt  and  clay  fill. 

g.  Impact  from  a  barge. 

h.  Shear  rupture  of  cross  wall  webs. 
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32.  The  variability  in  N ,  interlock  strength,  is  shown  by  the  histogram 
in  Figure  6,  which  is  based  on  a  Large  number  of  rests  by  Kay  (1975).  It  is 
obvious  that  the  histogram  is  skewed.  The  truncation  at  the  upper  end  is 
associated  with  web  failures  rather  than  interlock  failures,  and  it  is  likely 
that  a  greater  degree  of  symmetry  wouLd  have  been  produced  had  there  been  no 
web  failures.  There  is  some  justification  for  fitting  a  normal  distribution 
to  an  interlock  histogram  (Kay  1975).  Kay  has  indicated  that  a  normal 
distribution  with  a  coefficient  of  variation  (ratio  of  standard  deviation  to 
mean)  of  0.11  wilL  provide  a  good  representation  of  the  interlock  strength. 
Therefore,  a  deterministic  analysis  of  strength,  as  shown  alove,  results  in 
the  mean  value  of  the  interlock  strength;  then  a  value  of  0.11  can  be  taken 
for  the  coefficient  of  variation. 


INTERLOCK  STRENGTH  -  kips  per  inch 

Figure  6.  Measured  distribution  of  interlock  strength  (Kay  1975) 
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using  the  technique  described  in  Part  II,  is  presented  in  Figure  8.  This 
figure  also  presents  reliability  as  a  function  of  time,  t,  using  Lacroix's 
assumption  that  =  0.0025  in./yr  (Lacroix,  Esrig,  and  Luscher  1970). 

Shear  Rupture  in  Cross  Wall 

39.  A  cross  wall  may  become  so  corroded  that  it  will  rupture  in  shear. 
Since  the  hoop  tension  is  transmitted  to  the  Y-pile  at  the  junctions  between 
the  cylindrical  walls  and  the  cross  walls,  horizontal  tension  will  be  devel¬ 
oped  in  the  cross  walls  (Figures  9a  and  9b).  The  combination  of  a  cross  wall 
and  the  segments  of  the  corresponding  cylindrical  walls  act  as  an  I-beam  to 


TIME,  YEARS 

20  40  60  80  100  120 


6,  in.  (CORROSION  RATE  =  0.0025  in/yr.) 

Figure  8.  Estimated  condition  of  interlocking  system,  O'Brien  Lock  and  Dam 
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resist  the  shear  developed  by  the  lateral  pressures  on  both  sides  of  the  land 
wall  (Figure  9c).  Consequently,  an  element  of  the  cross  wall  near  the  junc¬ 
tion  with  the  Y-pile  is  subjected  to  shear,  S,  and  tension,  N'  (Figure  9d). 

By  the  Mohr's  Circle,  the  maximum  shear  (the  "load")  in  the  cross  wall  is: 

S  = _/s2  +  j  N'2  (19) 

max  v  4 

N'  is  the  tension  transmitted  to  the  cross  wall  and  is  expressed  as: 

N1  -  2N  cos  a  (20) 

where  N  is  the  hoop  tension  and  a  is  as  shown  in  Figure  9a. 

40.  This  shear  is  not  carried  by  the  cross  wall  alone;  however,  due  to  the 
lack  of  a  reliable  theory  to  estimate  how  much  of  the  shear  is  resisted  by  the 
fill,  it  is  assumed  that  all  of  it  is  carried  by  the  cross  wall.  For  the 
I-beams  with  heavy  flanges,  the  shearing  stress  is  nearly  constant  on  a  cross 
section  of  the  web  of  an  I-beam: 


S  =  P/b  (21) 

where  b  is  shown  in  Figure  9c,  and  the  definition  of  P  follows. 

41.  The  land  wall  of  the  lock  carries  greater  shear  than  the  river  wall 
because  of  the  pressure  of  the  backfill.  The  largest  total  shear,  P,  occurs 
at  the  bottom  of  the  lock  chamber  and  is  determined  from: 

P  =  L  JdHpdx  (22) 

where  p  is  the  net  lateral  pressure  acting  on  the  crosswall.  At  the  bottom  of 
the  chamber,  the  following  expression  can  be  substituted  for  p: 

p  =  Ka  y  (H  -  s)  +  (K  y'  +  Yw)  (s  -  x),  d  <  x  <  h  (23) 
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and  the  reliability  of  this  structual  element  is: 


R(e)  =  P(MS  >  0). 


(27) 


Figure  10  illustrates  this  reliability  function,  developed  for  the  shear 
strength  of  cross  walls. 


Tension  Rupture  of  Tie  Rod 


44.  A  guide  wall  is  formed  by  interlocking  Z-shaped  sheet  piles  driven 
into  clay  and  gLacial  till.  The  wall  is  partially  supported  by  a  system 
consisting  of  wales  and  tie  rods  (Figures  11  and  12).  The  tie  rods  are 


TIME,  YEARS 


0.00  Q05  O.IO  0.15  0.20  0  25  0.30 

£,  in  (CORROSION  RATE  =  0.0025  in./yr.) 

Figure  10.  Estimated  condition  of  cross  wall  shear  resistance, 

O'Brien  Lock  and  Dam 
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several  yards  apart.  One  bay  of  the  wall  between  successive  tie  rods  is 
considered  in  Equation  28.  The  tension,  P,  in  a  tie  rod  must  balance  the 
forces  on  one  bay  of  the  wall.  This  force  has  been  derived  from  Equation  25 
in  Kearney  (1986)  and  is  defined  as: 

P  =  \  [Kay(H2  -  s2)  +  (Kay*  +  Yw)s2  -  Ywh2  -  KY'd2]  (28) 


The  remainder  of  the  variables  are  defined  as  follows: 

L  =  distance  between  two  adjacent  rods  (84  in.) 

Ka  =  coefficient  of  lateral  earth  pressure  (variable) 
H  =  elevation  of  top  of  wall  (552  in.) 
s  =  elevation  of  backfill  saturated  water  (variable) 
h  =  elevation  of  water  in  river  (variable) 
a  =  elevation  of  wale  (444  in.) 


K  is  as  described  in  Equation  29. 


-K  y(H^~  s'*)  -  (K  y’  +  Y  )s2  +  y  h  2  +  KY’d"*  +  3K  ya(H2  -  s2) 

A  A  W  W  A 


+  3(K  y  ♦  Y  )as  -  3y  ah  -  3Ky 

A  w  w 


2  -  3Ky ' ad2  =  0 


(29) 


The  rod's  strength  is  calculated  from: 


„  .  _  u  (D  -  e) 

P  =  A  o  =  - -  o 

u  u  4  u 


(30) 


where  D  is  the  diameter  of  the  rod  (D  =  3  in,),  e  is  the  amount  of  corrosion, 

and  a  is  the  ultimate  tensile  strength  of  the  steel.  The  reliability 

u 

function  is  defined  as: 


R(e)  =  P(MS  >  0)  =  P(P  -  P  >  0)  (31) 

u 

Table  4  lists  the  random  variables  affecting  the  reliability  function  of  the 
tie  rods.  Figure  13  illustrates  the  reliability  function.  Although  the  tie 
rods  corroded  less  than  the  sheet  piles  themselves  (Kearney  1986),  the  lower 
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stability  is  considered.  A  deterministic  formulation  has  been  developed  for 
this  problem  (Kearney  1986).  The  loading  is  defined  as: 

ML  =  P(a  -  x)  -  \  (K  y'  +  y  *  K  y)(s  -  x)^  -  K  y(H  -  x)^ 

L  b  A  w  A  b  A 

+  t  Y  (h  -  x)"^  +  -jr  Ky  '  ( d  -  x)^  (32) 

b  w  b 

where  x  is  an  arbitrary  elevation,  P  is  computed  from  Eq  28,  and  K  is  computed 
from  Eq  29.  When  solved  for  the  point  of  maximum  moment,  x  is  as  follows: 

x  =  2[~(K.y'  +  Y  “  K.y)s  -  KyH  -  Y  h  +  Ky ' d ] / ( K  -  K  ) y '  (33) 

max  A  w  A  A  w  A 

On  the  other  hand,  the  resisting  moment  is  determined  by: 


M 


R 


o 

y 


l  - 


a  w 

y 


25E(t  -  e)' 


(34) 


where : 

0^  =  yield  point  of  the  steel 

w  =18  in.  (Figure  11) 

E  =  modulus  of  elasticity 

t  =  initial  thickness  (0.5  in.) 

I  =  moment  of  inertia  (18.37  cu  in.) 

Ke  =  the  effect  of  corrosion  on  the  moment  of  inertia  (55  cu  in.) 
C  =  5.5  in.  (Figure  11) 


Defining  the  reliability  function  as: 


R(e)  =  P(Md  -  M  >0)  (35) 

R  max 

and  using  Table  5  for  the  values  of  random  variables  yields  Figure  14,  which 
shows  the  effect  of  corrosion  on  the  stability  of  the  sheet  piles. 
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would  be  needed  to  understand  the  nature  of  this  problem  and  its  conse¬ 
quences.  These  findings  are  consistent  with  the  deterministic  service  Life 
assessment  in  Kearney  (1986).  However,  the  probabilistic  findings  are 
believed  to  be  more  useful  since  they  estimate  gradual  deterioration  as  an 
explicit  function  of  corrosion  and  time. 
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Figure  15.  Typical  cross  section  of  sheet  piling  bulkhead 


Load  Analysis 


50.  Sheet  pile  bulkheads  are  one  of  the  most  popular  types  of  water-soil 
supporting  systems,  although  corrosion  of  the  sheet  piles  seriously  affects 
the  bulkheads'  bending  moment  capacity.  However,  the  key  question  is  how  the 
actual  soil  load  distribution  is  developed.  Figure  16  illustrates  the  lateral 
earth  pressure  diagram  proposed  for  the  flexible  anchored  bulkheads  in  clean 
sand,  and  corresponds  to  the  results  of  the  model  test  with  flexible  anchored 
bulkheads  conducted  at  Princeton  University  (Tschebotariof f  1955).  This 
experiment  showed  that  full  restraint  of  the  lower  portion  of  the  bulkhead  was 
effective  when  the  ratio  D/H  (depth  [buried]  to  height)  equaled  0.43,  where 


TIE- BACK  POINT 


(a)  Assumed  deflected  shape  of 
sheet  pile  bulkhead 


(b)  Distributed  lateral  earth-and- 

water  pressure  on  fictitious  beam 


Figure  17.  Sheet  pile  analysis 


the  tensile  strength  of  the  sand  layer  that  has  been  saturated  by  the  capil¬ 
lary  action  above  the  water  level.  Because  of  these  considerations  and  the 
large  degree  of  uncertainty  associated  with  the  backfill  material  (it  is  not  a 
clean  sand),  in  these  calculations  has  a  mean  value  of  0.4  and  a  standard 
deviation  of  0.125. 

51.  The  one  common  factor  for  all  sheet  pile  bulkhead  structures  is  that 
the  maximum  bending  moment  occurs  in  the  underwater  zone.  Considering  this 
factor  and  measuring  down  from  the  water  level  (on  the  land  side),  the  loca¬ 
tion  of  zero  shear,  z,  is  determined  from: 


A  +  4  y  a2  +  y  az  -  4  K  yh2  -  K  yhz  -  4  K.y'z2  =  0 
p2w  w  2A  A  2  A 

or 

4  k.y'z2  +  (K  Yh  -  y  a)z  +  4  K.yh2  -  4  Y  -  A  =0 
2  A  A  w  2  A  2  w  p 


(38) 

(39) 
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at  that  time  is  t  -  e,  where  tQ  is  the  initial  thickness.  The  mean  section 
modulus  of  the  pile  per  unit  width  is  then: 

S  =  S  -  Ke  (42) 

o 

where  S  is  the  section  modulus  of  the  sheet  pile,  S0  is  the  initial  value  of 
the  section  modulus,  and  the  reduction  factor,  K,  is  a  constant  which  can  be 
determined  exactly  from  the  shape  of  the  cross  section.  K  can  also  be 
calculated  approximately  from  the  properties  of  the  cross  section.  The  result 
i  s : 


K  s  28.6  sq  in.  (43) 

Consequently,  the  mean  bending  capacity,  Mn,  of  this  sheet  pile  is  determined 
f  rom: 

M  =  o  (10.725  -  28. 6e)  (44) 

R  u 

Due  to  the  approximation  involved  in  determining  K,  a  variation  of  0.2  is 
chosen  for  M^,  while  a  mean  value  of  24,000  psi  with  a  standard  deviation  of 
480  psi  is  used  for  o^.  The  mean  value  of  is  calculated  from  Equation  41. 

Reliability  of  Sheet  Pile  Bulkhead 

53.  The  reliability  of  the  sheet  pile  bulkheads  is  investigated  by  consid¬ 
ering  a  margin  of  safety,  MS,  which  is  defined  as: 

MS  =  Mr  -  Ml  (45) 

Failure  occurs  when  this  margin  of  safety  becomes  negative.  MS  is  a  random 
variable  because  both  M^  and  M  are  random  variables.  The  reliability  function 
is  defined  as: 


R(e)  =  P(MS>0] 


(46) 
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PART  V:  PORT  STRUCTURES  IN  ISRAEL 


Corrosion  of  Steel  Sheet  Piles 


55.  BusLov  (1983)  has  developed  and  proposed  an  analysis  of  the  remaining 
service  life  of  sheet  piles,  applying  his  approach  to  the  port  structures  of 
Israel  that  were  inspected  for  corrosion  between  1976  and  1981.  Unlike  con¬ 
ventional  design  practice,  which  recommends  the  use  of  a  corrosion  allowance 
based  only  on  structural  considerations,  he  offers  the  addition  of  a  detailed 
analysis  of  the  actual  character  of  corrosion  to  the  design  process.  In  his 
approach,  the  sheet  pile  section  is  chosen  according  to  the  maximum  bending 
load  and  resistance;  however,  this  section  must  contain  some  corrosion  allow¬ 
ance.  In  the  area  of  active  corrosion  (mostly  in  the  tidal  and  splash  zone), 
the  existing  service  life  is  determined  by  considering  three  factors:  rate  of 
corrosion,  actual  bending  moments,  and  flange  thickness,  which  had  already 
been  determined  based  on  the  value  of  the  maximum  Dending  moment.  To  illus¬ 
trate  his  technique,  he  considers  a  corrosion  allowance  of  2.25  mm  (0.09  in.) 
based  on  a  30-year  life  span,  or  in  general,  corrosion  of  0.075  mm/yr.  This 
general  corrosion,  plus  the  maximum  bending  moment,  were  used  during  design  to 
choose  the  sheet  pile  section.  He  then  considers  the  upper  portions  of  the 
section  where  the  corrosion  is  rather  high  (Figure  19)  and  determines  the 
potential  service  life  for  several  corrosion  rates.  Figure  20  contains  the 
result  of  his  calculations  for  the  sheet  pile  bulkhead  at  Kishon  (Larssen  IVn 
Section;  t  =  0.58  in.). 


Probabilistic  Analysi s 

56.  The  current  research  considers  the  same  sheet  pile  bulkhead  in  Kishon, 
which  consists  of  a  concrete  deck  on  concrete  piles,  fronting  a  sheet  pile 
retaining  wall  (Buslov  1983).  Two  cases  are  considered:  corrosion  of  the 
sheet  pile  under  the  water  (corrosion  is  relatively  low  but  moments  are  rather 
large)  and  corrosion  in  the  upper  sections  of  the  sheet  piles  (corrosion  is 
usually  high  but  bending  moments  are  rather  small).  The  sheet  pile  cross 
section  is  a  Larssen  IVn.  The  corrosion  rate  is  considered  to  be  at  least  0.4 
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PART  VI:  SUMMARY  AND  RECOMMENDATIONS 


59.  The  reliability  function  of  an  existing  structure  serves  as  an 
index  of  its  structural  condition.  It  explicitly  reflects  changes  in  resis¬ 
tance  that  have  accumulated  during  service  as  well  as  changes  in  loading 
imposed  by  operational  conditions. 

60.  This  report  has  documented  the  development  of  a  straightforward, 
generally  applicable,  three-step  procedure  for  evaluating  structural  relia¬ 
bility.  The  margin  of  safety  is  assumed  to  be  normally  distributed,  and  point 
estimates  are  used  to  calculate  the  moments;  these  are  standard  probabilistic 
methods.  Its  computat ional  details  resemble  familiar  design  procedures  and 
are  thus  attractive  to  practicing  civil  engineers. 

61.  The  procedure  has  been  applied  to  three  different  types  of  sheet 
pile  structures  representative  of  those  maintained  by  the  Corps.  Emphasis  has 
been  placed  on  using  the  available  and  familiar  formulation  but  in  a  proba¬ 
bilistic  manner.  The  results  indicate  that  rehabilitative  cost  avoidances  may 
be  achieved  if  the  procedure  is  used  in  lieu  of  traditional  deterministic 
practices . 

62.  The  three-step  procedure  developed  herein  is  recommended  for  use  in 
evaluating  the  safety  of  existing  sheet  pile  and  other  metal  structures  oper¬ 
ated  and  maintained  by  the  Corps  of  Engineers.  The  reliability  of  structures 
as  calculated  by  this  procedure  should  be  used  as  a  component  of  the  condition 
index  of  the  information  system  under  development  to  manage  REMR  activities  in 
Corps  field  offices. 

63.  The  procedure  could  be  adapted  to  the  particular  characteristics  of 
concrete  structures  to  establish  a  consistent  evaluative  condition  index  for 
these  features. 
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Table  3 


Random  Variables  Affecting  Reliability  of 
Interlocking  System 


Variable 

Mean 

Standard  Deviation 

ka 

0.45 

0.125 

s ,  in . 

468 

20 

Y ,  Ib/cu  in. 

0.0723 

0.08  yw  =  0.0029 

Table  4 


Random  Variables  Affecting  Reliability  of  Tie  Rods 


Variable 

Mean 

Standard  Deviation 

ka 

0.45 

0.129 

h ,  in . 

444 

20 

Y ,  1 b/cu  in . 

0.0723 

0.0029 

s ,  in . 

520 

20 

R 

1.4 

0.333 

o  psi 

60,000 

6000 

Table  5 

Random  Variables  Affecting 

Rel iabi 1 ity 

of 

Sheet  Piles  Against 

Buckl ing 

Variable 

Mean 

Standard  Deviation 

E,  psi 

30,000,000 

3,000,000 

Oy,  psi 

36,000 

3,600 

Y  .  1  b / c u  ft 

0.0723 

0.0029 

ka 

0.45 

0.129 

R 

1.4 

0.333 

s ,  in . 

520 

20 

h,  in. 

444 

20 
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APPENDIX  A:  MICROCOMPUTER  PROGRAM 
FOR  ESTIMATING  THE  RELIABILITY 
OF  EXISTING  STRUCTURES 


This  appendix  provides  a  microcomputer  program  for  estimating  the  relia¬ 
bility  of  existing  structures.  The  program  is  written  in  BASIC  for  an  IBM  PC 
or  compatible  computer.  Translation  onto  other  hardware  or  into  other  soft¬ 
ware  is  straightforward.  The  main  program  in  lines  10  through  500  implements 
the  general  procedure.  The  subroutine  in  lines  1000  through  1130  specifically 
addresses  the  safety  margin  for  the  tutorial  miter  girder  example.  The  data 
statements  in  lines  2000  through  2020  contain  the  input  values  for  the  ex¬ 
ample.  Other  problems  can  be  solved  by  replacing  the  subroutine  and  the  data 
statements.  The  execution  of  the  program  for  the  tutorial  example  follows  the 
program  listing. 


A1 


1  t.nX>  51  lE'tRGIJ  I  I  NL  RE  fURNii  .  =  M:-j  F  OR  M  [  TEE  GIRDER  E 
1 0 1 0  A -41 . 3 2 
1  ,v20  hK-  '0.81 
1  30  D— ft  /’ .  b8 

1  •40  IN-  '■  1  7,v0 
1  /i'O  l.  740 

i  1  'fti  >  T  -ft  _  . 68 

i  1  i  '  '  7- 
•  'tin  F  V  -■  X  ■.  1  • 

|  ..'90  Ul  -  A  \  2  ) 

1  lO‘»  F  i  i  W .  AR  *  v  L  /  2  *  C  O  +  7  ) 

l!  to  F  B--W*  <-  I~  -  A 1  /  2 /_  I N*  a  L  *L  /  4  L»A»D'JMf-A'  '2-fl*fi) 
1120  v  -i  <  t  A  FB 
t  I  30  RETURN 

2000  DATA  "Miter  Girder  E.-ampl  e"  .  2,  "MS" 

20 1 o  DA T  A  "ty  psi "  , 39600 , 44o0 
..i  )20  DATA  "w  1  b  /  i  n  "  ,  1  Ouu ,  200 


A3 


tank  t: 


APPENDIX  B:  VARIABLE  LIST 


f 


a 


f 


bl 


a 


I 


w,  L,  Q , t 


MS 


x 


X 

X 


o 


X 


MS 


o 


MS 


ms 


R 


N 


u 


axial 

flexural  stress 

area  of  girder 

moment  of  inertia  of  girder 

Figure  3 

margin  of  safety 
arbitrary  elevation 

generic  input  random  variable  (Eq  4  through  10) 

mean  of  X 

standard  deviation  of  X 

mean  of  margin  of  safety 

standard  deviation  of  margin  of  safety 

margin  of  safety  calculated  using  mean  value  of  each 
variable 

reliability  (Eq  11) 

radius  of  circular  cell  (Eq  15)  (O'Brien:  316  in.) 
interlock  strength 
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elevation  of  the  top  of  the  fill  (O'Brien:  504  in.) 
(Eq  16) 

elevation  of  top  of  wall  (O'Brien:  552  in.)  (Eq  28) 
total  unit  weight  of  the  cell  fill  (variable) 
submerged  unit  weight  of  the  cell  fill  =  y  “  Y 

w 

unit  weight  of  water  (0.0361  Ib/cu  in.) 
coefficient  of  horizontal  earth  pressure  (variable) 
the  amount  of  corrosion 

time 

shearing  stress 
maximum  shearing  stress 

tension  transmitted  to  the  cross  wall 

Figure  9a 

Shear  force  in  a  cross  wall  bottom  of  a  chamber 
probability  (Eq  18) 

tension  in  a  tie  rod  (Eq  28) 

Figure  9c 

distance  top  of  bulk  head  to  point  of  tie  back 
attachment  (Eq  37) 

net  lateral  pressure  acting  on  the  cross  wall 
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the  effect  of  corrosion  on  moment  of  inertia 
O'Brien:  55  cu  in.) 
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Figure  11  (O'Brien:  5.5  in.) 

AP 

anchor  pull 

H,  b,  h,  L,  a 

Figure  17  (Eq  37) 

f" 

effect  of  wall  friction  on  reduction  of  active 

earth  pressure 

f ' 

uncertainty  in  relative  importance 

pressure  above  anchor,  and  tensile 

saturated  by  capillary  action  above 

of  passive  earth 

strength  of  sand 

water  level. 

z 

location  of  zero  shear 

M 

maximum  moment 

S 

section  modulus 

So 

initiaL  value  of  section  modulus 

K 

reduction  factor 

S(  e ) 

reflects  the  effects  of  corrosion  on  section  modulus 
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APPENDIX  C:  GLOSSARY  OF 
STATISTICAL  TERMS 


CUMULATIVE  DISTRIBUTION  FUNCTION:  The  cumulative  distribution  function  is 
denoted  F(x),  and  has  the  properties: 

( 1 )  Lim  F(x)  =  1 ; 

X-i-co 

( 2 )  Lim  F( x)  =  0 
x ■*“ 

(3)  F(x)  is  a  nondecreasing  function. 

MEAN:  The  sample  mean  or  average  of  a  set  of  n  measurements  x  ,  x2,...,xn  is 
the  sum  of  these  measurements  divided  by  n.  The  mean  is  denoted  by 
x,  which  is  expressed  operationally 
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NORMAL  DISTRIBUTION:  A  normal  distribution  has  a  bell-shaped  density 

_  (*-y)2 
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with  a  mean  =  p  and  a  standard  deviation  =  o. 

The  probability  of  the  interval  extending 

(1)  one  standard  deviation  each  side  of  the  mean: 

P[p  -o<x<p+oj  =  0.683 

(2)  two  standard  deviations  on  each  side  of  the  mean: 

P[p  -  2o  <  x  <  u  +  2oj  =  0.954 

(3)  three  standard  deviations  on  each  side  of  the  mean: 

P[u  -  3o  <  x  <  p  +  3o]  =  0.997 
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VARIANCE:  The  sample  variance,  s,  is  the  variation  of  individual  data  points 

about  the  mean.  The  dispersion  of  a  set  of  n  measurements  x  ,  x,,..., 

x  is  defined  as 
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. (x . -  x) 2 
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